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ABSTRACT 

The 10th recorded outburst of the recurrent eclipsing nova U Sco was observed simultaneously in 
X-ray, UV, and optical by XMM-Newton on days 22.9 and 34.9 after outburst. Two full passages of 
the companion in front of the nova ejecta were observed, witnessing the reformation of the accretion 
disk. 

On day 22.9, we observed smooth eclipses in UV and optical but deep dips in the X-ray light curve 
which disappeared by day 34.9, then yielding clean eclipses in all bands. X-ray dips can be caused 
by clumpy absorbing material that intersects the line of sight while moving along highly elliptical 
trajectories. Cold material from the companion could explain the absence of dips in UV and optical 
light. The disappearance of X-ray dips before day 34.9 implies significant progress in the formation 
of the disk. 

The X-ray spectra contain photosphcric continuum emission plus strong emission lines, but no clear 
absorption lines. Both continuum and emission lines in the X-ray spectra indicate a temperature 
increase from day 22.9 to day 34.9. We find clear evidence in the spectra and light curves for Thompson 
scattering of the photospheric emission from the white dwarf. Photospheric absorption lines can be 
smeared out during scattering in a plasma of fast electrons. We also find spectral signatures of resonant 
line scattering that lead to the observation of the strong emission lines. Their dominance could be a 
general phenomenon in high-inclination systems such as Cal 87. 
Subject headings: novae, cataclysmic variables - stars: individual (U Sco) 



1. INTRODUCTION 

A nova outburst is a thermonuclear explosion on the 
surface of a white dwarf star. The energy is produced by 
fusion of hydrogen to helium and is released via complex 
radiation transport processes through the optically thick 
ejecta, spanning a broad wavelength range. As the den- 
sity of the ejecta decreases with time, the peak of the ra- 
diation output continuously shifts from long wavelength 
optical emission to shorter wavelengths. The decreasing 
density allows us to see deeper into the outflow where 
the temperature and thus the energy of the emitted radi- 
ation is higher. Nova explosions occur in binary systems 
because the required amount of hydrogen-rich burning 
material can only be obtained by accretion from a sec- 
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ond star that orbits the white dwarf at a close enough 
dist ance to allow mas s transfer via an accretion disk, see 
e.g., iRobinsonl (|1976f) . Once all the previously accreted 
hydrogen is consumed, the nova will turn off, gradually 
returning to the original condition of the system. A nova 
outburst does not significantly change the system param- 
eters, and mass transfer can set in again, leading to new 
accumulation of hydrogen-rich material. Typically, the 
accumulation of sufficient hydrogen for another outburst 
is expected to take ~ 10 4 — 10 5 years, except for a hand- 
ful of recurrent novae, in which this cycle occurs fast 
enough that more th an one outbur st has been observed 
in a human lifetime (Schaefer 2010f). 

During the initial blast, the accretio n disk is destroyed, 
which was recently shown for U Sco by iDrake fc Orlandol 
(2010) utilizing hydrodynamic models. One question of 
interest is how quickly it reforms. This is not well known, 
and opportunities of direct observations of the reforma- 
tion process are rare. Significant theoretical efforts have 
been undertaken to understand the formation of accre- 
tion disks. An illustration of the expected formation pro- 
cess via Roche Lobe overflow is shown in fig. 1 in a re- 
view paper by iVerbuntl (|1982l ) . Initially, the gas streams 
with supersonic velocity along a highly eccentric orbit. 
As the stream hits the incoming gases, the flow down- 
stream changes into a circular orbit, which is expected 
to occur on time scales longer than radiative cooling and 
the dynamical or orbital timescale. In most cases, vis- 
cosity plays an important role for the details, but this 
fundamental picture remains the same. More important 
should be the presence of intensive high-energy radia- 
tion during the times of the XMM-Newton observations 
of U Sco. The recurrent outbursts in the eclipsing sys- 
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tern U Sco are a unique laboratory to observe the earliest 
stages of the disk formation. 

U Sco is a recurrent nova with now 10 recorded out- 
bursts since 1863. The latest outburst was discovered by 
B. Harris on 201 Jan 28.4385 (=HJD 2455224.94346) 
at V = 7.8 mag ()Schaefer et al.ll2010al ). iSchaefer et"aH 
(|2010bft extrapolated the light curve to earlier times, in- 
dicating that the peak magnitude was V = 7.5 on 2010 
Jan 28.1 (=JD 2455224.605), and that the time of out- 
burst was 



to = 2010Jan27.8 (= JD 2455224.305). 



(1) 



This time is used as the reference time to in this arti- 
cle. The unde rlying sy s tem w as identified as an eclips- 
ing system by|Schacfer (199Q|) with an inclination angle 
of > 80° (|Thoroughgood et al.ll200"l . Based on 29 ac- 
curate eclipse times from 2001-2009, we find a best fit 
linear ephemeris in units of days to be 



H J Drain = 2451234.5387 + E * 1.23054695, 



(2) 



where E is an integer that counts t he cycles . Th e 
binary separation is 6.5 ± 0.4 R (|Schaeferl 119901) . 
and the component masse s have been determined by 
iThoroughgood et al.l (|2001h to be > 1.31 M Q for the 
white dwarf and 0.88 ± 0.17M Q for the secondary. The 
radius of the companion is 2.1 ± 0.2 R Q , indicating that 
it is evolved (|Schaeferlll990T ). The high mass of the white 
dwarf, close to the Chandrasekhar mass limit, is con- 
sistent with the short recurrence time scale and led to 
spec ulations that USco m ay be a supernova la progen- 
itor dStarrficld et al.lll98S ). Ho wever, recent abundance 
measurements bv lMasonl (1 2 Oil indicate that the under- 
lying white dwarf may of the ONeMg kind which does 
not contain enough nuclear binding energy for a SN la 
exposion. Consequently, even if the Chandrasekhar mass 
limit is reached, the white dwarf may turn into a neutron 
star via core collapse. 
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Fig. 1.— Swift/XKT (top) and UVOT (bottom) light curves of 
U Sco. The shaded regions mark the duration of continuous, deeper 
Chandra (ii) and XMM-Newton (t2 and 43) observations. 

Swift X-ray a nd UV monitoring sta rted within a day 
after discovery (|Schlegel et al.l[2010aL but there was no 
X-ray detection until 2010 Jan 31 (day 4). The Swift X- 
ray and UV light curves shown in Fig. [T] will be discussed 
in more detail by Pagnotta et al. (in preparation) . Until 
~day 8, the X-ray light curve was essen tially constant 
with weak, hard emission above 3 keV (Schlcgel et al. 



I2010al) . On ~day 12, the Swift count rate was found 
to be 100 times higher , exhibiting a supersoft spectrum 
(|Schlegel et al.l l2010bl) that resembles the class of su- 
persoft sources (SSS, Kahabka fc van den Heuvel 119971 : 
Ivan den Heuvel et al.1 11992P) . At the same time, until 
~day 12, the UV flux dramatically decreased, supporting 
the picture of a shrinking photospheric radius, accompa- 
nied by a continuous increase in photospheric tempera- 
ture. Between days 12 and 25, no changes in X-ray and 
UV brightness were seen, while after day 25, another 
increase in X-rays with a simultaneous decrease in UV 
emission level occurred. From day 14.25 - 17.54 after dis- 
covery, two well covered, and one p artly covered, eclipses 
were seen in the UV light curve bv lOsborneetaLl pOlOh 
who also reported that, in X-rays, no clear, sharp, deep 
or total eclipses were seen. The X-ray light curve was 
highly variable with typically lower flux around the time 
of eclipse. Around 2010 Mar 8 (day t + ~ 40 days), 
the brightness of U Sco experi enced a sharp drop in the 
X-ray/UV/Optical/IR bands (|Schaefer et al.ll2010cD . in- 
dicating that nuclear burning had turned off. 

Reports of deeper and continuous Chandra and XMM- 
Newton observations taken on 2010 February 14 and 19 
(t% = t + 18.7 and = tg + 22.9 days) , respectivel y were 
given bv lOrio et all (|2010l) and iNess et ail (f2010h . The 
times of these observations plus a second XMM-Newton 
observation taken on t% — to + 34.9 days are marked 
by gray shaded areas in Fig. [TJ bracketing the respec- 
tive start and stop times. High-resolution X-ray grating 
spectra consist of photospheric continuum emission plus 
broad emissio n lines that varie d greatly between the two 
observations (|Ness et al.H2010D . yielding a systematic in- 
crease in flux in lines of higher ionization stage. 

2. X-RAY/UV OBSERVATIONS 

Two 16-hour XMM-Newton observations were carried 
out with start times 2010 February 19.65 {t% = to + 
22.9 days) and 2010 March 3.7 (t 3 = t + 34.9 days). Si- 
multaneous X-ray light curves (0.1-1 keV), X-ray spectra 
(5-38 A range, 0.01 A resolution), and ultraviolet (UV) 
light curves and spectra (2200-3600 A range, 15 A reso- 
lution) were extracted from the European Photon Imag- 
ing Camera (EPIC), the Reflection Grating Spectrome- 
ter (RGS), and the Optical Monitor (OM), respectively. 
For day 22.9, 27 UV grism spectra (approximately 2000 
seconds exposure time) were taken and were used to ex- 
tract UV light curves by integrating over three different 
wavelength intervals: the full range, 2200-3600 A, and 
the ranges 2800-3600 A and 2200-2800 A. In addition, 
we extracted an uncalibrated zero-order light curve from 
the non-dispersed photons which represent the brightness 
over a broad wavelength band with central wavelength 
- 4000 A. Since the zero order is not flux calibrated, 
count rates are only presented as an indicator for relative 
brightness variations. For day 34.9, only band-integrated 
fluxes are available but in high time resolution of the OM 
fast mode. 

Separate X-ray light curves were extracted in a soft 
(0.1-0.5 keV) and a hard (0.5-1 keV) band, S and H, and 
were used to compute a spectral hardness light curve 
HR= (H - S)/(H + S). 

3. LIGHT CURVES 
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TABLE 1 
Observation log 
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In Fig. [21 a time series of X-ray brightness, X-ray hard- 
ness, UV flux, and zero-order (optical blue) count rate 
are shown as a function of orbital phase, based on Eq. [5J 
The optical and UV light curves (bottom two panels, er- 
ror bars are smaller than the plot symbols) contain clean 
primary eclipses. The optical light curve is smoother 
than the UV light curve. The X-ray light curve (top 
panel) exhibits high-amplitude variations between two 
brightness levels (horizontal green and orange lines) . The 
hardness ratio shown in the second panel is less variable 
than the X-ray brightness. A broad, flat peak may be 
seen that is centered around phase 1.0, indicating that 
softer emission from the central source may be eclipsed. 
The event at phase 1.21 could be a flare-like brightening, 
as anomalies are seen in all bands plus in the spectral 
hardness light curve; for more discussion on the spectral 
evolution, see Sect. 14.1.31 Possibly, also the sharp in- 
crease in X-rays plus reduction in spectral hardness at 
phase 0.9 could indicate a flare, but no simultaneous 
signatures in the UV/optical or hardness light curves 
support that. Note that similar sudden brightness in- 
creases have been seen in other novae, e.g., V1494Aql 



([Drake et al.ll200l . 

After phase 1.25 (quadrature), the count rate stabi- 
lized at the high-flux state for the remaining ~ 20 ks 
(=5.5 hours), ~ 30% of the entire observation. This phe- 
nomenon is worth closer inspection because after phase 
1.25, the X-ray dark secondary moves behind the X-ray 
bright photosphere which can then no longer be eclipsed. 
To investigate whether this could be a coincidence, we 
searched for other X-ray observations taken outside of 
eclipse, but found no similar occurrence in other data. 
The large majority of Swift observations (see Fig.[l| were 
centered around eclipses, and no continuous light curve 
outside eclipse can be assembled from these data that 
could c onfirm or reject this conclusion. It is notewor- 
thy that lOsborne et al.l (|2010D report typically lower flux 
around the time of eclipse, based on Swift/XRT data cov- 
ering several eclipses. The Chandra light curve extracted 
from the observation on t± = to + 18.7 days was taken 
outside ecl i pse, a nd no such variations were reported by 
lOrio et al~l (|2010h . but this observation was quite short 
(15 ks) and much earlier in the SSS phase. Although co- 
incidence cannot be completely ruled out, we consider 
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Fig. 2. — Short-term evolution of X-ray and UV/optical brightness and X-ray hardness (second panel) after <2 = to + 22.9 days, with 
bandpasses given in vertical axes and legends. Time units are phase, based on Eq. [2] Phase 1.0 is when the companion is between the 
white dwarf primary and the observer, and the vertical dotted line indicates quadrature phase. After phase 1.25, the companion is located 
behind the primary. In the top horizontal axis, hours relative to phase 1.0 are given. Simultaneous Swift data are added with blue star 
symbols, after scaling by a factor 25. UV grism fluxes are plotted in the third panel in units of 10 — 15 erg cm -2 s^ 1 A -1 . 



the fact that one third of the entire observation shows 
no dips as support for the conclusion that the variations 
are related to the binary orbit and therefore, that the 
occulted plasma resides inside the binary orbit. 

Except for the two events at phases 0.9 and 1.21, we 
consider the patterns of variability as flux decrements 
(dips) with the undisturbed flux level being that observed 
after phase 1.25. An inhomogeneous brightness distri- 
bution that is eclipsed by the companion can securely 
be excluded, because during four minima the flux level 
drops to 50%. Each of four bright regions would have to 
contribute 50% to the total flux which is not possible. 

The similar flux levels during minima and maxima 
leads us to the conclusion that the same source is oc- 
culted multiple times by different occultcrs. Since dur- 
ing the minima the flux is not zero, the central source 
can not contribute more than ~ 50% to the total X- 
ray emission, while the rest must originate from further 
away, e.g., caused by scattering in the outer regions of 
the ejecta (see also Sects . I4TTT1 and 14" . 1 . 3[) . 

Within this picture, the totally different behavior in 



the UV and optical can be explained by occulting plasma 
that is transparent to UV, only leaving visible a single 
eclipse by the companion. For example, neutral hydro- 
gen is transparent to UV/optical light but opaque to soft 
X-rays. Alternatively, the X-ray emitting material that 
is likely more concentrated to around the white dwarf 
suffers sharper occultations by clumps while the larger 
optical emitting region will not have any significant frac- 
tion covered. 

During the second XMM-Newton observation on day 
34.9, the nova was already in the process of turning off, 
as is evident from a rapid decline i n X-ray and UV bright - 
ness observed by Swift (Fig.Q]and Schacf er et al.l l2010c). 
This can also be seen in the XMM-Newton light curves 
in Fig. [31 where the X-ray and UV light curves are shown 
in the top and bottom panels, respectively. Swift/XRT 
count rates, scaled by a factor 25 to compensate for dif- 
ferent sensitivity, are also shown. The X-ray decline rate 
is marked by a dotted line with text above it. After con- 
version for different sensitivities, this rate is consistent 
with the long-term decline rate in the Swift/XRT light 
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Fig. 3. — Evolution of X-ray and UV brightness (top and bot- 
tom) plus X-ray hardness (middle), taken on t$ = to + 34.9 days, as 
in Fig. [2] The dashed line in the top panel marks the X-ray decline 
rate of 6.2 cps day -1 and is consistent, after rescaling by the XRT- 
to-EPIC scaling factor of 25 (see caption of Fig.[2l, with the long- 
term decline rate in the Sm/i/XRT data (see Fig.Q. In the bottom 
panel, UV fluxes from XMM-Newton/ OM and Swift/UVOT are 
shown in units of 10 — 15 erg cm -2 s — 1 A -1 . In both light curves, a 
small step is seen at phase ~ 0.6 and the light curves are asym- 
metric, yielding a sharper egress than ingress. This could indicate 
an occultation of the central source by an accretion stream that is 
located ahead of the path of the companion. 

curve that is marked in Fig. [T] 

The Swift/UVOT fluxes are consistent with the OM 
fluxes (bottom panel). A small step at phase ~ 0.62 can 
not be explained by an eclipse by the companion, but it 
could be an eclipse by an accretion stream (see $5|). 

A clean eclipse can be seen in all light curves centered 
around phase 1. The hardness light curve shown in the 
middle panel does not change with the eclipse. Towards 
the end of the observation, the source becomes slightly 
softer, possibly related to the decline. 

The dips have disappeared from the X-ray light curve, 
and must therefore be understood as a transient phe- 
nomenon during the early SSS phase. This behavior ap- 
pears similar to RS Oph and other novae that have been 
monitored with Swift, showing an early variability phase. 
The XMM-Newton data of U Sco give valuable insights 
into the processes that may explain such an early vari- 
ability phase in more general terms. 

3.1. Light curve model 

As a qualitative test of the interpretation of occulta- 
tions of the central source by dense absorbers, we have 
constructed a geometrical model which consists of two 
sources of light and four occulters. On a 350 x 350 bright- 
ness pixel map, two sources of light are defined, a non- 
variable, spherical central source with a radius i? pr i ma ry 
plus additional constant emission that is not subject to 
occultations. Light originating from the central source 
is removed by the secondary star or by three clumps of 
gas that is optically thick to X-rays. All occulters are as- 
sumed to be spherical. In this way, a synthetic light curve 
can be generated by calculating for each phase bin the 
difference between the sources of light and the light orig- 
inating from those pixels located behind the absorbers 
as seen from Earth. 

The inclination angle and radii of the orbit and com- 
panion are well known and are fixed parameters. Un- 
known parameters are the sizes of the central source and 
the absorbers, the location of the absorbers, and their 



opacity (see next paragraph) . Values of these parameters 
were obtained by a combination of manual and automatic 
iteration with the aim to obtain good agreement between 
the synthetic and the observed light curves. Many dif- 
ferent configurations are possible such as a non-spherical 
geometry of the source or the absorbers as well as in- 
trinsic variability. Under these circumstances, we do not 
believe that beyond a qualitative feasibility test, reliable 
quantitative conclusions can be derived. 

The opacity of the absorbers is parameterized by the 
column density of neutral atomic hydrogen, N-r. De- 
pending on this value, a certain fraction of soft X-ray 
emission that originates from behind a given absorber 
can penetrate, yielding full transparency for Ah = 0. 
The fraction of light that can pass through the absorber, 
T(wi), depends on wavelength w;: T(wi) = e - N H><<?(wi) ; 
where o~(wi) is the cross section for a given wave- 
length bin Wi and was computed usi ng the bamabs tool 
which is part of the Pintof Ale package ( Kashva p fc Drakel 
2000). The fraction of absorbed light, integrated over 
the entire spectrum, sp(wi), is F = 1 — ^2 i {T{wi) x 
sp{w i )}/y d ,sp{w i ). 

In Fig. 21 a graphical illustration of a best- fit model 
to the X-ray light curve from day 22.9 is shown in the 
top. Below, the resulting light curve in comparison to 
the observed X-ray light curve and the residuals relative 
to the measurement errors (x) are shown. Except for 
the flare- like event at phase 1.2 and the sharp rise at 
phase 0.9, all residuals are consistent with the degree 
of variability after phase 1.25 and can be considered as 
not related to occultations. Further refinements of this 
model will not increase the precision of the conclusion. 

While the deg ree of reproduction is not great, xi-cd — 
4.7, we attribute the deficiencies to additional random 
intrinsic variability that can not be parameterized. The 
parameters of this model are, for the central source, a ra- 
dius of -Rprimary = 0.67 R© (17% of the Roche Lobe) and 
a brightness level of 16.6 counts per second. A power-law 
surface brightness distribution with a power law index a 
was assumed, but the index came out very small with 
a = —5.8 x 10 -5 . We experimented with limb darkening 
or -brightening but found no better reproduction of the 
data. The level of additional constant emission was 17.7 
counts per second. The nature of the additional non- 
variable emission is discussed in the context of the spec- 
tra in Sect. I4.1TT1 and [4.1.3l The absorbers are assumed 
to orbit at a radius of 3.55 R©, having a size of 0.73, 1.1, 
and 0.83 R©, respectively. For the column density in the 
absorbers, a lower limit of Ah > 1.4 x 10 21 cm~ 2 is de- 
termined. This is illustrated in Fig. [5] from which it can 
be deduced that Ah is unconstrained at high values. We 
conclude that the absorbers are completely opaque to the 
X-ray spectrum of U Sco, and the lower limit represents 
the value above which only the constant additional emis- 
sion remains. This can be compared to predictions from 
the expected density of the accretion disk material and 
the path length through it. The volu me of a flared disk 
rim is of the order of 10 34 cm 3 , e.g., iDrake &: Orlandol 
(IMol). For a mass loss rate of ~ 10 _7 M© yr _1 (e.g., 
lHachisu et al.l[2000[ ). and a build-up time of ~ 5 days, the 
accretion disk rim density is of order 3 x 10 14 cm~ 3 . For 
a path length through the gas of ~ 0.1 R©, the column 
density is of order 10 24 cm~ 2 . While this is higher than 
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Fig. 4. — Qualitative geometrical model to the X-ray light curve 
taken on day 22.9. The top panel gives a spatial representation 
of the model as seen from Earth (inclination angle given in top 
right legend) and is to scale with units given in the axis labels. 
The location and size of the central source is marked by the black 
filled circle, centered at the origin of the coordinate system and 
the assumed radius given in the top left legend. The central source 
is eclipsed by the companion (red circle) and by three co-rotating 
absorbers (blue, green, and orange circles). The sizes and posi- 
tions of the circles are to scale with the solid and dotted styles 
denoting the respective positions at the start and the end of the 
observation. Eclipses by any of these occulters lead to temporary 
decreases in total light that, in the absence of eclipses, is assumed 
to be constant. Additional emission is added that remains con- 
stant, regardless of binary phase. The resulting model light curve 
is shown in comparison to the data in the middle panel. The bot- 
tom panel shows the difference between data and model, relative 
to the measurement errors, \ =(data-model)/A. Except for two 
events at phases 0.9 and 1.2, the residuals between model and data 
are of the same order as the variations outside eclipse, past phase 
1.25. 



the number derived from the eclipse model, a fraction of 
all hydrogen can be ionized, requiring a higher column 
density for the same opacity, or the build-up time may be 
shorter. Moreover, the disk may not have fully reformed 
at the time of observation, and during the reformation 
process, the density should be lower. 

The vertical dotted line in the bottom two panels of 
Fig. [1] indicates phase 1.0, where the center of eclipse by 
the companion should be. A dip of the expected width 
is present, however, it is clearly shifted. For better re- 
production of the data, we had to include a phase shift 
parameter with a value of -0.026 in the model. Such an 
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Fig. 5. — Determination of column density of absorbers in light 
curve model. While the light model yields poor resemblance with 
observations for low values of TVjj, any high value of TVjj is consis- 
tent with the data. 

artificial shift could be avoided if a non spherical symme- 
try of the central X-ray source is allowed. For example, 
we tested a model with a crescent shaped source and were 
able to reproduce the second dip without a phase shift. 
There is, however, no physical explanation for a crescent 
shaped central source. 

The same procedure was applied to the 27 OM grism 
exposures of day 22.9. We constructed two separate 
UV light curves by integrating each grism spectrum over 
the two wavelength bands 2200-2800 A and 2800-3600 A. 
Combined with the light curve from the non-dispersed 
photons, we have three light curves representing different 
spectral colors, all containing clear eclipses. In Fig. [6] the 
results are illustrated, where the red circles represent the 
location and relative size of the companion during each 
of the 27 exposures and the shaded circles in the cen- 
ter represent the primary with the best-fit radii for each 
band, following the color scheme given in the top right 
legend. The best-fit radii are given in the top left leg- 
end. The short-wavelength UV band yields the smallest 
radius, while the longest-wavelength optical band (zero 
order) yields the largest radius, close to the inner La- 
grangian point. The significance of the light curve fit is 
illustrated in Fig. where a range of primary radii is 
plotted against changes in x 2 relative to the best fit for 
each wavelength band. The boundaries around the 1- 
(7 uncertainties are marked by the vertical colored lines, 
based on an increase of x 2 by 1 (horizontal black line). 
The best-fit light curve models are shown in comparison 
to the data (in normalized units) in the bottom panel of 
Fig. [6] While the models (points connected by lines, see 
bottom right legend) give no perfect reproduction of the 
data (shaded areas and gray line, see top left legend), 
the width of the respective eclipses are well reproduced, 
yielding reliable radius estimates. 

The trend of decreasing radius with decreasing wave- 
length appears like a temperature gradient. Hotter 
plasma, producing more short-wavelength emission, re- 
sides closer to the center than cooler plasma. 

4. SPECTRA 

The X-ray spectra from two types of spectrometers are 
discussed in tj4.ll High-resolution soft grating spectra 
and low-resolution hard CCD spectra are available for 
both observations. In addition, 27 UV grism spectra 
have been taken during the first observations that are 
discussed in £14.21 

4.1. The X-ray spectra 
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Fig. 6. — Geometrical light curve model applied to three 
UV/optical light curves extracted from the non-dispersed zero 
grism order (optical, orange) and two UV light curves (green and 
blue) , extracted from different wavelength ranges as denoted in the 
upper right legend. The central source was assumed spherical with 
a power-law brightness distribution and the resulting best-fit radii 
for each wavelength band given in the top left legend. The three 
emitting regions are illustrated as filled circles. The location and 
size of the companion is plotted as a red circle at the beginning 
of each of the 27 exposures. In the bottom panel, the three ob- 
served light curves (shades) and the best-fit models (points) are 
shown in normalized units. 



XMM-Newton carries two types of X-ray spectrome- 
ters that give information from different energy bands, 
yielding different sensitivity and spectral resolution. The 
EPIC consists of two types of CCD detectors, the MPS 
(Metal Oxide Semi-conduct or; I Turner et "all 120011 ) and 
the pn (|Striider et all 120011) . Two identical MOS cam- 
eras, MOS1 and MOS2, are in operation behind two 
separate mirrors, and they can be operated in different 
modes. They are sensitive between 0.1-10 keV with an ef- 
fective area ranging between 50 — 500 cm 2 and an energy 
resolution of roughly 50 eV. The pn has a lower resolu- 
tion of 100 eV but a higher sensitivity of 200 - 1300 cm 2 
from 0.1 — lOkeV because of the larger collecting area of 
an extra set of mirrors. 

The dispersive spec trometers RGS1 and RGS2 
(jden Herder et all 120011 ) share the light with the same 
mirrors as the MOS cameras. They are both sensitive 
between 6 — 38 A (0.3 — 2keV) with an effective area be- 
tween 20 — 60 cm 2 and a wavelength resolution of 0.05 A. 
The dispersed photons are recorded by MOS type CCDs. 
The intrinsic energy resolution of the detectors is used to 
separate higher dispersion orders from the first disper- 
sion order spectrum. While the sensitivity of the RGS is 
much lower than the EPIC, the superior resolving power 
of AE/E ~ 100 - 600 compared to AE/E < 1 for EPIC 
allows individual line transitions to be identified. For ex- 
ample, only in the RGS spectrum can an absorption line 
spectrum be distinguished from an emission line spec- 
trum, provided the source is bright enough or the expo- 
sure time is sufficiently long for enough signal to noise. 
SSS spectra are extremely bright and are therefore ide- 
ally suited for studies with the RGS. The EPIC spectra 
are superior above ~ 2 keV and are use to study the hard 
Wien tail and to search for harder emission that may be 
related to shocks in the surrounding medium, or accre- 
tion. 

Separate spectra for different time intervals can be ex- 
tracted by application of filters on photon arrival times. 
We have extracted the RGS and EPIC spectra from the 
entire observations and a series of spectra from adjacent 
shorter time intervals in order to probe for spectral time 
evolution. Times were converted to orbital phase using 
Eq. 1 
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Fig. 7. — Estimate of uncertainties in the primary radius obtained 
from three UV/optical light curves, using the color scheme defined 
in Fig. [6] (see legend). 



4.1.1. High-resolution RGS spectra 

The merged RGS1+RGS2 spectra from both observa- 
tions are shown in photon flux units in Fig. |8] The spec- 
tra shown here have been produced with the SAS com- 
mand rgsf luxer. The spectral range of the RGS covers 
the peak of the absorbed photospheric emission from the 
optically thick ejecta. A combination of continuum emis- 
sion and strong emission lines can be identified. 

As a parameterization of the continuum component, 
we fitted a blackbody model to those spectral bins that 
are not dominated by strong emission lines. The result- 
ing models are marked by the yellow/blue curves in Fig. [8] 
for each observation (see caption) . The resulting param- 
eters effective temperature, interstellar neutral hydro- 
gen column density, and effective radius (derived from 
dilution factor and assumed distance 12kpc) are given 
in the left legends. The bolometric luminosity, derived 
from these values applying the Stefan-Boltzmann law, 
are also given but an unrealistically high value is derived 
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Fig. 8. — RGS spectra taken on day 22.9 (top) and 34.9 (bottom). 
Flux units are phcm -2 s -1 A -1 . A blackbody curve is included as 
blue/yellow line, where those spectral bins that were used for fit- 
ting are yellow. The best-fit parameter values are given in the left 
legends but the main purpose of the blackbody fitting is to show 
that the continuum component is of atmospheric origin. The pho- 
tosphcric temperature has greatly increased from day 22.9 to day 
34.9, evidenced by a clear shift of the Wien tail to shorter wave- 
lengths. Additional strong emission lines are marked with labels. 
Lines of higher ionization stages are seen on day 34.9, yielding 
additional evidence for a temperature increase. 

for day 22.9, exceeding the Eddington luminosity by two 
orders of magnitude. This is a well-k nown deficiency 
of bla ckbody fitting to SSS spectra (e.g.. iKrautter et all 
1996), and we therefore caution that these values can 
not be trusted to reflect the true physical conditions, 
and they are only given for reproducibility. Appropriate 
models have to account for radiation transport in non- 
local thermal equilibrium (NLTE), plus the e xpansion of 
the e.jecta has to be taken into account (e.g.. |]NTessll2010t 
Ivan Rossum fc Nesil2010l ). These are currently unsolved 
problems, and their discussion is beyond the scope of this 
paper. Nevertheless, the reproduction of the continuum 
component by a blackbody model is sufficiently good to 
allow the conclusion that the continuum originates from 
the photosphere around the white dwarf. The absence of 
any absorption lines could indicate that we are not di- 
rectly seeing the atmospheric emission. Gray scattering 
in a plasma of fast electrons (Thompson scattering) could 
explain this observation, as narrow atmospheric features 
would be smeared out. 

While the photospheric temperature can not be 
accurately determined from blackbody fits, the location 
of the Wien tail clearly indicates that the tempera- 
ture must be well above 10 5 K. From this lower limit 
of the temperature we can conclude that the radial 
extent must be significantly smaller than the radius 
of the Roche Lobe (4.1 R©). Otherwise, the bolo- 
metric luminosity would be four orders of magnitude 
above the Eddington limit which is physically unrealistic. 

The emission lines originate from highly ionized nitro- 
gen, oxygen and carbon. For Nvn and Nvi, the H-likc 
and He-like ls-np series can be identified up to S (ls-5p). 
In the second observation, additional Ne ix and Ne x lines 
can be identified. 

In order to escape, the emission lines must originate in 
a reasonably optically thin plasma, which can reside at 



any distance, far above the photosphere. On the other 
hand, the strongest lines are seen at wavelengths where 
the continuum is strongest, and they are thus likely pho- 
toexcited, requiring that they come from hot plasma re- 
siding close to the ionizing source. 

This type of X-ray spectrum, in both observations, 
is not typical of novae such as V4743 Sgr, RS Oph, or 
V2491 Cyg in which pro minent deep photospheric ab - 
sorption lines were seen bv lNess et afl (|200ai20Q7ll20ll . 
It resembles more that of the ecl ipsing supersoft source 
Cal87 which lEbisawa et al.l (|2010l) explain as a combi- 
nation of Thompson scattered continuum plus emission 
lines from resonant line scattering as viewed at high in- 
clination angle (see also iGreiner etaLll200l . The pho- 
tosphere around the white dwarf cannot be seen directly 
as it is blocked by the accretion disk. A similar situ- 
ation may be occurring in U Sco and perhaps even in 
all high-inclination systems as wa s already sugge stive 
from UV observations pre sented by lLa Dousl (|1991l ) and 
iVitello fc Shlosmanl (|1993|) . Thompson scattering in the 
electron-rich environment would preserve the spectral 
color, but a non-zero velocity distribution of electrons 
would smear out photospheric absorption lines which can 
explain the absence of any clear indication of absorption 
lines in the RGS spectra of U Sco. 

In addition to Thompson scattering, resonant line scat- 
tering may occur, selectively amplifying or diminishing 
resonance emission lines, depending on the geometry of 
the scattering plasma. For example, the geometry of 
an accretion disk allows more photons from resonance 
transitions to be scattered into the line of sight than 
out of the line of sight, and resonance lines can appear 
amplified when compared to forbidden lines. This is 
best seen in the He-like Nvi triplet at 28-29 A, where 
the resonance line is stronger than the intercombina- 
tion line. For comparison, in the high-resolution X-ray 
spectra of low-density stellar coronal plasma, the for- 
bidden line is usually almost equally as strong as the 
reso nance line and st ronger than the intercombination 
line (|Ness et al.l 120011 ). While the low ratio of forbidden 
to intercombination can be interpreted as high densities 
(|Gabriel fc Jord an 196i|), the same low ratio can also 
be caused by strong UV radiation fields, if the emitting 
plasma resides close to the photosphere. Inspection of 
dip and non-dip spectra presented in Sect . 14 . L~3l confirms 
this explanation. 

4.1.2. EPIC spectra 

EPIC spectra yield the best information about the 
high-energy part of the spectrum, above 2keV, where 
the RGS is not sensit ive. Swift detected h ard emission 
early in the evolution (jSchlegel et al.ll2~010al ). and it is of 
interest to use the deeper XMM-Newton EPIC data to 
check whether it was still present on days 22.9 and 34.9. 

In Fig. IH1 the photon image of MOS1 taken on day 
22.9 is shown with all events in the left panel and only 
those with registered energies above 2keV in the right 
panel. While in the right panel, a clear source seems to 
be present, this can not be interpreted as the detection of 
hard emission. The blue circle marks the radius within 
which significant pile up is present, which was determined 
using the SAS tool epatplot. All photons in the right 
panel are within this circle, indicating that they are con- 
taminated by pile up. We encountered the same problem 
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Fig. 9.— EPIC/MOS1 observation taken on day 22.9. In the 
left panel, the pixel positions of all events between 0.1-10keV are 
shown relative to the peak of the point spread function. In the 
right panel, only photons between 2 and 10 keV are shown. A 
concentration of counts towards the source position can be seen, 
however, the inner 180 pixels (marked by the blue circle around 
the center) are affected by pile up (see text). 




Fig. 10. — EPIC/MOS2 timing mode observation taken on day 
22.9. In this mode, all photons along the detector y axis (RAWY) 
are accumulated into a single number of counts for each detector 
x-pixel (=RAWX) coordinate. The source spectrum was extracted 
from between RAWX values 290-310 (marked by gray shaded area), 
encircling 74% of all counts. Plotted is the number of counts per 
pixel, integrated over three different energy ranges as indicated in 
the legend. Above 2keV, no significant detection is found. The 
energy range 1.25-1.5 keV contains the H-like and He-like Mg lines. 

in the pn observation and in all EPIC observations of the 
second observation. 

More sensitive to faint hard emission are the MOS2 
data in the first observation which was operated in tim- 
ing mode (see Table [JJ), yielding a shorter readout time. 
We determined, with epatplot, that pile up was not a 
problem. In Fig. [TOl the brightness profiles as a func- 
tion of the detector coordinates (RAWX) are shown for 
three different energy ranges as indicated in the legend. 
The full-range profile clearly indicates the location of the 
source in RAWX coordinates, and above 2keV, noth- 
ing is seen at the same position. In the narrow range 
1.25-1.5 keV, containing the H-like and He-like Mg lines, 
a clear detection is found, indicating that some harder 
emission above the SSS component is present. A source 
spectrum can be extracted from a range in RAWX of 
290-310 pixels, which is marked by the gray shaded area 
in Fig. [TUJ No other sources are located along the y axis 
that could contaminate the source region. For our further 
considerations, we therefore only use the MOS2 spectrum 
for the first observation while for the second observation, 
we use the most sensitive pn data, after correction for pile 
up. 

The EPIC spectra for days 22.9 and 34.9 are shown in 
Fig. [TTJ The energies of some line transitions that have 
been seen in other novae are marked with vertical lines 
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Fig. 11.— EPIC spectra taken on days 22.9 (top, MOS2 detector) 
and 34.9 (bottom, pn detector). The energies of strong H-like and 
He-like transitions are marked, and the red lines indicate blackbody 
models with the same parameters as used in Fig. El F° r day 22.9, 
a model with an additional optically thin collisional component 
(APEC) is plotted with a blue thick line. 

with labels. 

The SSS continuum appears featureless at EPIC en- 
ergy resolution. Blackbody models with the same pa- 
rameters as used in Fig. [8] are overplotted with the red 
curves. Obviously, the EPIC spectra are not suitable 
for studies of the soft part below 2 keV, even if they are 
not piled up. However, in the first observation, there is 
significant excess emission above the blackbody curve at 
energies Z 0.7 keV, indicating the presence of an addi- 
tional plasma component. This excess can be modeled 
with an optically thin thermal APEC model (Smit h et al.l 
120011 ) as indicated by the blue thick curve in Fig. [TTJ The 
APEC model component assumes collisional equilibrium 
for ionization and excitation processes which is observ- 
able as bremsstrahlung continuum plus emission lines. 
Optically thin emis sion has been seen in various no- 
vae be f ore, see, e.g.. ILTovd et al.l (l!992Tk iKrautter et al.l 
(fl996tk iMukai fc Ishidal (|200lD ; (Orio et all (j200lH and 
can be caused by shocks with pre-existing surround- 
ing material, e.g., from previous outbursts, or within 
inhomogeneous ejecta. The electron temperature ob- 
tained from fitting an APEC model plus a blackbody 
component with fixed parameters to the spectrum is 
kT = 0.29 ± 0.03 keV. This temperature corresponds to 
an average random electron velocity of 10 9 cms -1 when 
assuming a Maxwellian velocity distribution, which is ap- 
proximately twice the expansion velocity. The H-like and 
He-like Ne and Mg ls-2p lines are well reproduced by the 
APEC model component. The Ovm line at ~ 0.66 keV 
(18.97 A) is not reproduced by the combined bb+apec 
model and is more likely photoionized by the SSS con- 
tinuum. In the second observation the collisional com- 
ponent has faded significantly. 

Above 2 keV, we found no significant continuum emis- 
sion nor Fe fluorescent lines and lines of FeXXV at 6.4 
and 6.7 keV, respectively. A more detailed discussion of 
the evolution of the hard emission will be instructive in 
the context of all X-ray observations, including the early 
Swift and Suzaku data. 

4.1.3. Spectral Evolution 
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A) Long term evolution 

A comparison of the spectra from the two observations 
yields the spectral evolution on time scales of days to 
weeks. A tr end of inc re asing temperature has been 
reported by INess et all ([2010D based on increasing 
strengths in high ionization lines from day 18.7 to day 
22.9. A gradual increase of the effective temperature 
is expected as the photospheric radius continues to 
shrink during the SSS phase. From the spectra shown in 
Figs. 151 and [XT] this trend can be seen to have continued 
through day 34.9, evidenced by a clear shift of the Wien 
tail of the continuum to higher shorter wavelength and 
the stronger emission lines belonging to high-ionization 
stages (Fig. [8]). Only the high-ionization He/H-like Mg 
lines are stronger on day 22.9 than on day 34.9 (Fig. [IB"]) . 
Since they are not formed on top of any continuum 
emission, they can not be photoexcited. The more likely 
origin of the Mg and Ne lines is a collisional plasma, 
and this interpretation is supported by the APEC model. 

B) Short term evolution 

To study the short term evolution, we have extracted 
multiple RGS and EPIC spectra covering adjacent 
time intervals. In Figs. [T^] and [T21 the spectra and 
pn light curves, rotated by 90 degrees, are shown in 
the top and right panels, respectively. The central 
parts are brightness maps with increasing orbital phase 
from top to bottom and increasing wavelength/energy 
from left to right for the RGS and EPIC data, respec- 
tively. The key to the color scheme is given in the top 
right. In addition, the hardness light curves shown in the 
second panels of Figs. [5] and [3] are used for the discussion. 

To illustrate the effects of contamination of the high- 
energy part of the EPIC spectra by solar flare activity, 
field-integrated flare particle background light curves are 
shown with light blue color in Fig. 1131 To determine 
the particle background, a full-field light curve was ex- 
tracted above lOkeV. For better comparison, the parti- 
cle background light curve is amplified by a factor 20. 
Above 2 keV, no significant emission can be identified in 
the EPIC spectra of either observation. Around phase 
1.3 in the first observation and phase 0.85 in the second, 
some harder emission was present that could be contam- 
ination because it coincides with elevated levels of solar 
flare particle background emission. 

We have also extracted separate RGS spectra for 
times during dips and peaks for day 22.9 and in- 
and outside of eclipse for day 34.9. A combined 
dip spectrum has been extracted from the time in- 
tervals (in spacecraft units 10 8 s after 1998-01-01) 
3.8298385 - 3.8298715, 3.8298925 - 3.8299995, and 
3.8300645 — 3.8301275. A peak spectrum is integrated 
over the five time intervals 3.8298166 - 3.8298362, 
3.8298726 - 3.8298878, 3.8300018 - 3.8300620, 
3.8301286 - 3.8301858, and 3.8302500 - 3.8304502. 
The inside eclipse spectrum was extracted from the 
time interval 3.8404789 - 3.8407309 and the outside 
eclipse spectrum from 3.8401449 - 3.8403669. During 
extraction, the spectra were binned to 0.1-A bins, larger 
than the default 0.01 A for better signal to noise. The 
differences between the resulting spectra are shown in 
Fig. 031 The shaded areas around the difference spectra 



represent the uncertainty ranges derived via error prop- 
agation from the measurement errors in the individual 
spectra. If the times of low flux are interpreted as 
occultations of emission originating from within the 
binary orbit, the difference spectra are representative of 
the plasma residing in the inner orbit while the low-flux 
spectra represent only outer regions. 

Bl) Flare-like events 

Two significant brightness increases occurred around 
phases 0.9 and 1.2 in the first observation which in the 
following are referred to as flares. In the top panel of 
Fig. [T21 it can be seen in the contour levels that both 
events were accompanied by increased emission in the 
Nvi lines which originate from lower ionization plasma 
than the other emission lines. During the flare at phase 
1.2, some extra emission in the Nvi intercombination 
line can be seen. At the same time, the Nvil Lya line 
at 24.8 A continues to be fainter during the dip while the 
higher Ly series lines of N VII seem to peak around the 
same time as the Nvi lines. 

The contemporaneous appearance of the lower- 
temperature line with both flare-like events is unexpected 
because flares are usually high-temperature events in- 
volving some heating mechanism. For example, in solar 
and stellar coronae, flare events are usually accompa- 
nied by the appearance of high-ionization e mission lines 
which are produced directly in the flare fe.g. JGiidel et all 
12001 . 

When assuming photospheric events as origin of the 
flares and that the Nvi line emission has responded to 
these changes, then a time delay of ~ 0.01 in phase, 
corresponding to ~ 17 minutes after the peak of each 
flare, can be used to derive the maximum distance 
between the plasma that emits the excess in Nvi and 
the photosphere. The two black horizontal dashed lines 
in Fig. [12] illustrate the time delay, and it corresponds 
to a light travel distance of ~ 2 AU. This is outside of 
the binary orbit of 6.5 R© where cooler plasma resides. 
Hotter plasma further inside may have been blocked 
during the times of the flares, as both occurred during 
dips. 

B2) Spectral changes on day 22.9 (dipping) 

Closer inspection of the spectral changes with the dips 
in the first observation gives insights into the geometry of 
the emitting gas at this time. The brightness map in the 
top panel of Fig. [12] illustrates that during the dips, the 
continuum and the emission lines are both fainter, but 
neither component is completely occulted. Both com- 
ponents must therefore originate from both, the inner 
orbit and from further outside. Thompson scattering of 
the continuum component has already been discussed in 
Sect. 14.1.11 to resolve the conflict between high temper- 
ature from the Wien tail and the required small size of 
the central source plus the absence of absorption lines. 

A comparison of the emission lines during and outside 
dips on day 22.9 can be studied in Fig. [14] The difference 
spectrum represents the plasma within the inner orbit, 
revealing the spectral signatures of resonance line scat- 
tering. While the strong ls-2p resonance lines of, e.g., 
Nvn (24.78 A) and Nvi (28.74 A) originate only from 
the outer regions, lines belonging to transitions involv- 
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Fig. 12. — Illustration of time evolution of RGS spectra taken on days 22.9 (top) and 34.9 (bottom). The central part in each plot is a 
brightness map composed of 63 and 54 separate spectra (1000s exposure time), respectively, with colors representing increasing flux levels 
from light green, yellow, red, purple, dark blue, and light blue. The corresponding flux for each color can be determined using the vertical 
legend bar in the top right. In the top panel, two isocontour levels are included in white. Increasing wavelengths are plotted from left 
to right, and increasing phase from top to bottom; phases have been calculated from photon arrival times based on Eq. [2] The dashed 
horizontal lines running across the brightness map indicate selected times for which the corresponding spectra are shown with colored 
shades in the top panel. In addition, the average spectra from Fig. [8] are included. In the right, the pn light curves from the top panels of 
Figs.[2]and[3]are shown, rotated by 90° clockwise. The dips in the top panel and the eclipse in the bottom panel are seen in both continuum 
and emission lines in the brightness map. Flare-like events at phases 1.21 and 0.9 on day 22.9 yield an increase in the Nvi lines at ~ 29 A. 
The two black lines illustrate a time lag between flares and the changes in Nvi that correspond to a light travel distance of 2AU. 

ing higher principal quantum numbers are stronger in the difference spectrum, e.g., the Nvi He- like series lines 
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Fig. 13. — Illustration of the spectral evolution in the .2-2 .7 keV 
band in the same format as explained in the legend to Fig. 1121 based 
on 63 and 62 EPIC spectra, extracted from adjacent time intervals 
in the observations, starting day 22.9 (top, MOS2 detector) and 
34.9 (bottom, pn detector), respectively. The white lines in the 
brightness maps are isocontours. The energies of strong H-like and 
He- like transitions are marked in the top parts (see also Fig. Ill I) , 
In the right, the same pn light curve from the top panel of Fig. [2] is 
shown, rotated by 90° clockwise. In addition, the field-integrated 
light curves extracted from above lOkeV are shown to illustrate 
times at which solar flare activity may compromise the high-energy 
part of the spectra. 

He/3, 7, and 8. The respective oscillator strengths are 
/ = 0.17, 0.05, and 0.02, which are low compared to the 
oscillator strength of the Hea line of / = 0.71. Line pho- 
tons with high values of / can be absorbed and re-emitted 
in a different direction. In a non-spherical plasma, strong 
resonance lines can be amplified or reduced, compared to 
lines of lower oscillator strengths. The fact that we are 
seeing the high-/ lines from the inner regions reduced is 
consistent with the geometry of an accretion disk. 

Forbidden transitions can not be absorbed, and the 
contributions from the inner region to the ratio of the 
Nvi intercombination and forbidden lines at 29.1 and 
29.54 A is unaltered by resonance scattering. Owing to 
the proximity to the central source, the low f/i ratio re- 
ported in Sect. I4.1TT1 has to be interpreted as photoexci- 
tation rather than high density. 

In the top panel of FigfTBl the evolution of the Wien 
tail can be studied that is a qualitative indicator for the 
photospheric temperature. During dips, the continuum 
emission has gone down at all wavelength, and no shift 
of the Wien tail can be seen that would indicate that the 
continuum component from the outside regions has a 
different temperature. This supports the interpretation 
of Thompson scattering. The sudden change in the 
variability pattern at phase 1.25 is accompanied by a 
sudden increase in brightness above 0.6 keV, extending 
to ~ 1.5 keV, including the Mgxi ls-2p transition at 
1.3 keV. While solar flare activity could be held respon- 



sible for the brighter emission between phases 1.2-1.35, 
the spectra taken after phase 1.35 are not contaminated 
and are still harder. This could be interpreted in the 
general context of increasing temperature on longer 
time scales as discussed above. Also, a slight increase 
in the high-ionization O vm line is suggestive in the top 
panel of Fig. but it may also be part of the general 
increase in brightness. 

B3) Spectral changes on day 34.9 (eclipse) 

The spectral evolution during the second observation 
is displayed in the bottom panels of Figs . IT21 and [TBI The 
continuum and the emission line components experience 
similar changes during the eclipse as during the dips on 
day 22.9. The total emission decreases by ~ 50% during 
eclipse, indicative of significant scattering processes. The 
emission lines decrease by roughly the same amount as 
the continuum. 

This is illustrated with the difference spectra shown 
with blue color in Fig. [Mj The N VI and N vn lines are 
hardly present, and therefore, all the plasma emitting the 
observed N lines resides outside of the binary orbit while 
some of the flux in the hotter O vm line also originates 
from within the orbit. In order to explain the changes 
from day 22.9 to 34.9, the line emission that was first 
seen to come from the inner regions (Nvi He/3 and in- 
tercombination line) has to be moved to further outside 
to escape eclipses by the companion on day 34.9. While 
for the Nvi He-like intercombination line one could ar- 
gue that the entire N VI triplet is too faint for studies of 
changes, the line complex around 25 A is strong enough 
that residual Nvi He/3 would be detectable. Also, the 
higher Ly series lines of N VII are not seen in the differ- 
ence spectrum. We interpret the absence also of high Ly 
series lines in the difference spectrum as evidence for an 
increase in the optical depth of the scattering plasma, 
which could be part of the reformation process of the 
accretion disk. A bit puzzling appears the presence of 
O VIII Lya in the difference spectrum while Ly/3 is not 
present. 

Except for the Ovin Lya line, the reduction in 
brightness during eclipse can therefore exclusively be at- 
tributed to the continuum. The shape of the contin- 
uum in the difference spectrum in Fig. [14] is the shape 
of the continuum component in the inner regions. It has 
roughly the same shape as the out-of-eclipse spectrum, 
supporting the interpretation of achromatic Thompson 
scattering. 

Close inspection of the RGS and EPIC brightness maps 
in the bottom panels of Figs. [T2l and [TBI reveals that the 
eclipse progresses in a slightly non-uniform way in wave- 
length/energy. In the Ray leigh- Jeans tail, between 26- 
28 A, the continuum seems to go through a wider eclipse 
towards longer wavelength (Fig. [T"2"]) . In the Wien tail, 
the eclipse seems to be narrower towards higher energies 
which is indicated by the contours in Fig. [TBJ This could 
indicate a temperature gradient, however, the hardness 
light curve in the middle panel of Fig.[B]does not indicate 
a significant temperature change with the eclipse. Only 
the softening after the eclipse appears noteworthy and 
may have to be attributed to the cooling while nuclear 
burning is turning off. On the other hand, the hard- 
ness light curve contains the emission lines and is a less 
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Fig. 14. — Difference spectra between dip- and peak spectra for day 22.9 (gray) and in- and out of eclipse spectra for day 34.9 (light 
blue), representing the spectra of the missing emission during dip and eclipse, respectively, thus from the inner-binary region. The indices 
Lya, /3, and 7 indicate H-like Lyman series lines with ls-2p to ls-4p transitions, and Hea to S are the respective He-like series lines with 
transitions ls-2p to l-5p. The notation Hea, i, and / behind the Nvi labels indicate He-like ls-2p transitions with 1 S- 1 P (resonance line), 
(intercombination line), and 1 S- 3 S (forbidden line), respectively. The Nvi intercombination line (i at 29. 1 A) and the He-like series 
lines p to <5 originate from the inner orbit on day 22.9 and are thus photoexcited by the nearby photospheric continuum emission. The Nvi 
Hea line is reduced in the inner-orbit spectrum because of resonant line scattering out of the line of sight. The Ovm resonance line at 
18.97 A on day 34.9 reaches us more directly from the inner regions while the N lines are reduced by the same amount as the surrounding 
continuum. 



sensitive indicator for the photospheric temperature. 

4.2. UV spectra 

In the first of the two XMM-Newton observations, 
starting day 22.9, the UV grism of the OM was em- 
ployed, taking 27 consecutive spectra (see Tabled]). In 
Fig. 1151 these spectra are shown in the same format as 
in Figs. [12] and [[3] In the range 2200-3600 A, the spec- 
trum shows many features which are difficult to identify. 
Wavelengths of a few known lines are marked. A weak 
feature coincides with the position of the 2800 A Mgn 
resonance doublet but, given the presence of many sim- 
ilar features in the spectrum as well as the uncertainty 
of the absolute wavelength scale (which depends on the 
accurate centroiding of the zero order), its identification 
is uncertain. 

The wavelength range of the OM grisms overlaps with 
the long-wavelength range of the International Ultravi- 
olet Explorer (IUE). Five spectra were taken with the 
IUE during the outburst in 1979: on June 28, June 30, 
July 2, July 4, and July 10, corresponding to 5, 7 , 8, 
10, and 16 days after outburst (Will iams et al.lll98lD . In 
the top panel of Fig. [15] the IUE spectra from 8 and 16 
days after the 1979 outburst are overplotted in orange 
and blue (see right legend). The strongest lines in these 
spectra are, apart from the Mgn doublet, Hen at 2512, 
2734 and 3204 A, Neiv at 2421 A and [Mgv] at 2784 and 
2928 A. 

The eclipse minimum can be seen to coincide with 
reductions in flux at all wavelengths in the brightness 
image, but some spectral ranges, e.g., ~ 2350 A, yield 
narrower eclipses than others, e.g., ~ 3200 A. The two 
dashed black lines in the image draw attention to the pos- 
sibility of broader eclipses towards longer wavelengths as 
found in g3l]for the two wavelength bands 2200-2800 A 



and 2800-3600 A. If these bands are interpreted as spec- 
tral colors, the conclusion would be that hotter plasma 
resides further inside. However, the concept of spectral 
color as temperature indicator is only valid for contin- 
uum sources while emission lines can arise at any wave- 
length, independently of their formation temperature. 
While emission lines of higher ionization stages are pro- 
duced closer to the central ionizing source, they do not 
necessarily arise at shorter wavelengths. It can be seen 
that some features, e.g., at 2650 A, have indeed narrower 
eclipse profiles than others, e.g., at 3125 A. Also, not all 
features have their eclipse minimum at the same phase, 
suggesting that, if they are real lines, they arise from 
different emission regions in an inhomogeneous plasma. 

5. DISCUSSION 

The different behavior in the simultaneous XMM- 
Newton X-ray and UV light curves starting on day 22.9 
was unexpected with a clear eclipse only in UV and 
optical but multiple dips in the X-rays. Some similarities 
with the class of X-ray dippers (Walter et al.lll981l 119821) 
may be h elpful in discuss in g pos sible scenarios. For 
example, Di az Trigo et al.l (|2009l ) presented in their 
figure 5 the XMM-Newton X-ray and UV light curves 
of the high-inclination X-ray binary XB 1254-690 and 
found clear eclipses in the V band while in the X-ray 
band, deep dips were seen at phase 0.8. While multiple 
dips are generally not seen in the X-ray dippers, EXO 
0748-676 containe d dips at phases 0.6 , 0.7, 0.8, 0.9 and 
0.1 (see figure 9 in lParmar et al.l[l986| ). While a number 
of similarities are thus seen, USco is still a completely 
different system, and any parallels to X-ray dippers have 
to be interpreted with care. 

In X-ray dippers, significant changes in spectral 
X-ray hardness have oftentimes been observed, while 
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Fig. 15. — Illustration of the spectral evolution in the same format as Fig. 1121 for 27 OM U grism observations, starting on day 22.9. 
In the top panel, two IUE spectra taken 8 and 16 days after the 1979 outburst of USco are included in orange and blue as indicated in 
the right legend. In the right panel, the X-ray light curve and the OM light curve are shown in black and blue, respectively. The two 
diverging black dashed lines in the central panel are drawn to guide the eye, indicating a possible trend of broader eclipse towards longer 
wavelengths. The light blue dashed line indicates the time of the flare at phase 1.2. 



in U Sco, the hardness remained the same during dips. 
An achro matic ste e p dec rease in flux has also been 
seen by iSala et all (|2008l) in the nova V5116Sgr. In 
contrast to U Sco and other novae, the central object in 
X-ray dippers is an accreting neutron star, emitting a 
much harder continuum. Harder X-rays can partially 
penetrate material of large column densities, and only 
the soft tail of the X-ray spectrum is blocked. The 
continuum spectrum of U Sco is much softer, and any 
material of significant column density above ~ 10 21 cm~ 2 
will block all emission from behind, regardless of energy. 
Total disappearance of all X-ray light during dips is 
only avoided because of Thompson scattering of the 
central emission into a much larger region rather than 
partial transparency of the absorbers. Our parameter 
studies of iVn in the light curve model support this 
interpretation (see Fig. [5]). An example of achromatic 
dipping at low energies as a consequence of high- Ah 
absorbers in relation to the shape of the sou rce spectrum 
is Cyg X-2 (|Balucinska-Church et all 120111) . While the 
continuum is much harder than in U Sco, the opacity of 
the absorbers is ~ 10 23 cm -2 , leading to full blockage 
of all light in the 0.5-2 keV band. The hardness ratio 
from 0.5-1.0/1.0-2.0 band is therefore constant during 
the dips, while only emission from the 4.5-8.0 keV band 
can partially penetrate. In U Sco, not enough emission 
at energies above 1 keV is present to constrain Ah of the 
absorbers in the way this is possible for Cyg X-2. 

A question of interest is the reformation of the 
accretion disk that was destro yed during the nova 
outburst ([Drake fc Orlandol 120101 ) . The transition from 



dipping on day 22.9 to clean eclipse on day 34.9 may 
be part of the reformation process, and the concept of 
intersecting opaque clumps as part of the reformation 
process deserves closer consideration. The theoretical 
pos sibility of cl umps in the ejecta has been discussed 
by iDiaz et all (|2010f ). however, being immersed in a 
high-radiation environment, the ejecta must be highly 
ionized, and as such, they do not absorb soft X-rays via 
the photoelectric effect. Thompson scattering in clumpy 
material can also be ruled out as cause for the dips as 
that would also lead to the same dips in the UV/optical 
light curves. Clumpy ejecta can therefore not cause 
such deep dips in the X-ray light curve. A possible 
source of opaque clumps consisting of cooler material 
is the companion star, if accretion has set in again. 
Cooler plasma is opaque to soft X-ray while at the same 
time being transparent to light with wavelengths longer 
than 912 A (13.6 eV), the ionization energy of atomic 
hydrogen, which would explain why no dips are seen in 
the simultaneous UV and optical light curves. Absorb- 
ing atomic material from the companion is thus strong 
support for the interpretation that accretion has set in 
again. This does not require the full reestablishment 
of the accretion disk, a clumpy accretion stream can 
also be imagined. The absence of dips after phase 1.25 
suggests that no clumpy material is on the far side as 
seen from the companion which is an argument in favor 
of an accretion stream, residing between the companion 
and the central white dwarf. 

In Fig. [TBI we show a theoretically calculated initial 
trajectory of a reforming disk based on equations 1 and 
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Fig. 16. — Trajectory of accretion stream at an early ti me of the 
reformation of the accretion disk. The calculation follows Flannervl 
( 1975) using the parameters of U Sco. The trajectory interacts with 
itself, forming an asymmetric ring at the beginning of disk refor- 
mation. The radius of the central source is drawn assuming the 
radius derived from the light curve model (see top panel of Fig. (3) , 
which is also consistent with the blackbody radius. The trajectory 
dives into the optically thick ejecta on the far side from the inner 
Lagrangian point (LI). This could explain the residual variability 
seen in the bottom panel of Fig. [J] Being inside the photosphere, 
this part of the accretion stream does not occult the photosphcric 
X-ray emission. No dips are therefore expected between phases 
1.33 and 0.75, while during other times, dips can occur, depending 
on the structure of the accretion stream. At a later time, the circu- 
lar ring that is the protodisk interacts with the incoming stream, 
leading to circularization that could lead to the disappearance of 
dips in X-rays as seen on day 34.9 (see Fig[3j|. 

I in IFlannervl(|T975l) . assuming the parameters of U Sco. 
During the early disk formation phase, an asymmetrical 
"ring" is formed. The phases that correspond to differ- 
ent viewing angles are marked. The hashed area in the 
center encircles the optically thick parts of the ejecta. 
The surface can be considered the source of X-ray 
emission, and no views inside are possible. The highly 
elliptical trajectory runs through the optically thick 
part of the ejecta, and any inhomogeneities inside the 
photosphere can not cause dips. An eccentric trajectory 
of infalling material can thus be the explanation for not 
seeing any dips after phase 1.25. 

A remarkabl y simi lar event has been observed by 
lin't Zand et al.l ([20111 ) in the accreting neutron star in 
2S 0918-549. After an X-ray burst, strong (87 ± 1% 
peak-to-peak amplitude) achromatic fluctuations were 
seen as a transient phenomenon. Since the outburst 
evolved much more rapidly, the X-ray light curve taken 
with the Proportional Counter Array (PCA) on the Rossi 
X-ray Timing Explorer (RXTE) may be interpreted as 
a fast-mo tion picture of the ey olution in USco (see fig- 
ure 1 in lin't Zand et al.l |20TT1). It was speculated by 
m 't Zand et al.l (|2011h that the fluctuations were due to 
Thompson scattering by inhomogeneities in a resettling 
accretion disk that was disrupted by the effects of super- 
Eddington fluxes, much like in USc o where the blast 
wave has destroyed the accretion disk ([Drake fc Orlandol 
2010). Furthermore, they argue that an expanding 



shell may be the necessary prerequisite for the fluctu- 
ations. Irregularities in the accretion stream could be 
caused by the continued expansion of the nova ejecta 
that has stopped between days 22.9 and 34.9. The stream 
could be disrupted by the ejecta, or the companion star 
may be distort e d in the outer layers. A difference to 
lin't Zand et al.l ([20111 ) is the interpretation of the prop- 
erties of the intervening material, yielding highly ion- 
ized plasma. Dips are caused by Thompson scattering 
out of the line of sight rather than photoelectric ab- 
sorption. This is consistent with their observation of in- 
creases above the decay trend of the burst which can be 
explained by sideways or backward scattering of clouds 
located behind the compact object. Such increases are 
not seen in U Sco. Furthermore, achromatic Thomp- 
son scattering would also leave its footprints in the UV 
and optical light curves, and photoelectric absorption by 
high- Ajj plasma in a low state of ionization appears more 
likely for U Sco. 

6. SUMMARY AND CONCLUSIONS 

From the simultaneous XMM-Newton X-ray, UV, and 
optical observations of the 10th recorded outburst of the 
recurrent, eclipsing nova U Sco, we have deduced new in- 
sights into time variability, spectral properties, and spa- 
tial structure: 

• Time variability: As opposed to expectations 
that eclipse mapping will reveal the photospheric 
radius of the white dwarf in different stages of 
the evolution, we have seen achromatic dipping 
in X-rays but clean eclipses in UV and optical. 
Dips were only seen before quadrature, phase 1.25. 
The later observation contains clean eclipses in all 
bands. During dips and eclipses, ~ 50% of X-ray 
emission remained, indicating that the observed X- 
ray emission is extended beyond the binary orbit. 

• Spectral Properties: The high-resolution X-ray 
spectra consist of blackbody-like continuum emis- 
sion plus strong emission lines. No absorption lines 
can be identified. From day 22.9 to day 34.9, an 
increase in temperature can be deduced from the 
shape of the continuum and from increases in emis- 
sion lines arising in ions of higher ionization stages. 

The eclipsing supersoft source Cal87 yields a 
qualitatively similar high-resolution X-ray spec- 
trum. Photospheric X-ray emission from the cen- 
tral source is blocked by the accretion disk and can 
only be seen after Thompson scattering in the ac- 
cretion disk corona. Strong emission lines come 
from resonant line scattering within the surround- 
ing optically thin medium. Only emission lines 
from transitions with low oscillator strengths are 
seen to originate directly from the inner regions. 

The UV grism spectra on day 22.9 contain emission 
lines that undergo eclipses, which indicates that 
they originate from within the inner binary orbit. 
The interpretation as an emission line spectrum is 
supported by a comparison to the 1979 outburst, 
during which the UV spectra were already optically 
thin 7 days after outburst. 

• Evolution of the Spatial Structure: We inter- 
pret the dipping as the result of occultations of the 
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central source by high-density absorbing gas that 
is aligned along the trajectory of a reforming ac- 
cretion stream. The gas consists of plasma of low 
ionization stage, originating from the companion 
star. The continued outflow disrupts the accre- 
tion stream into fragments that intersect the line 
of sight. After nuclear burning has ceased, the out- 
flow stops, no more fragmentation takes place, al- 
lowing the accretion stream to circularize and con- 
dense into a fiat accretion disk, yielding no more 
dips. 

One important argument in favor of recstablishment 
of accretion is the low degree of ionization of the 
intervening material which indicates that it originates 
from the companion. The most important supporting 
observations are the achromatic nature of dips and the 
lack of dips in the other wavelength bands. Ionized 
plasma would not cause any photoelectric absorption 
while low-ionization plasma is only opaque to soft X-rays 
but transparent to optical/UV emission, longward of 
912 A. While the dips could be caused by Thompson 
scattering in the electron-rich environment of the highly 
ionized nova ejecta, the same dips would have to be seen 
in UV and optical. 

The observed spectral signatures for combined Thomp- 
son scattering and resonant line scattering are: 

• The continuum component does not go to zero dur- 
ing either dips or eclipse. The location of the 
Wien tail requires the source to be small enough 
for complete occultations to avoid super Edding- 
ton luminosity. Achromatic Thompson scattering 
can increase the geometric size of the emission re- 
gion without changing the shape of the continuum. 
Photospheric absorption lines are smeared out in 
the fast-moving scattering plasma and are there- 
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